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Abstract. We generalize some earlier results on a Berezin-Toeplitz type of 
quantization on Hilbert spaces built over certain matrix domains. In the 
present, wider setting, the theory could be applied to systems possessing sev- 
eral kinematic and internal degrees of freedom. Our analysis leads to an iden- 
tification of those observables, in this general context, which admit a semi- 
classical limit and those for which no such limit exists. It turns out that the 
latter class of observables involve the internal degrees of freedom in an intrin- 
sic way. Mathematically, the theory, being a generalization of the standard 
Berezin-Toeplitz quantization, points the way to applying such a quantization 
technique to possibly non-commutative spaces, to the extent that points in 
phase space are now replaced by N X N matrices. 



1. Introduction 

Let n be a symplectic manifold, with symplectic form u>, and H, a subspace of 
L 2 (f2, dfj,), for some measure \x. For <j> £ C°°(fi), the (generalized) Toeplitz operator 
Tcf, with symbol <p is the operator on TL defined by 

(1.1) T^f = P(<f,f), feTL, 

where P : L 2 (il, dfi) — > TL is the orthogonal projection. It is easily seen that T$ is a 
bounded operator whenever is a bounded function, and ||T^||^_,.-^ < ||0||oo, the 
supremum norm of (f>. 

Suppose now that both the measure fi and the subspace TL are made to depend 
on an additional parameter h > (shortly to be interpreted as the Planck constant), 
in such a way that the associated Toeplitz operators on TLh satisfy, as h \ 0, 

(1-2) \\T^\\n^H h - Halloo, 

and 

(1-3) WT^T^-T^Wn^n^O, 

(1-4) HIP? 5 ,T^]-T^ } \\n h ^n h ^0 

(where {*,*} is the Poisson bracket with respect to u>), and, more generally, 

oo 

for some bilinear differential operators Cj : C°°(Q) x C°°(f2) — ► C°°(f2), with 
C ((f),ip) = 4>ip and C\{(j),il)) — Ci(tp,<f>) = -^{(f>,ip}. Here the last asymptotic 
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expansion means, more precisely, that 
(1-6) 

^W-E^SU) « „ =0(h N+1 ) a S h\0, W = 0,1,2,.... 

r. fl-h * '*-h 

3=0 

One then speaks of the Berezin-Toeplitz quantization. Indeed, it is well known that 
the recipe 

oc 
j=0 

then gives a star-product on f2, and (jl.3p . (jTT4j) just amount to its correct semiclas- 
sical limit. 

The simplest instance of the above situation is £1 — R 2 ™ ~ C", with the standard 
(Euclidean) symplectic structure, and 

(1.7) H h =L 2 h jn,dn h ) 

the Segal-Bargmann space of all holomorphic functions square-integrable with re- 
spect to the Gaussian measure dfj>h(z) := e~" z ^ / h (irh)~ n dz idz being the Lebesgue 
measure on C n ). As shown by Coburn [Cob , (JT7SJ) then holds with 

(1.8) (</>»= V -,d a 4> • T$. 

\ a \=3 

The resulting star-product coincides, essentially, with the familiar Moyal product. 

Other examples of Berezin-Toeplitz quantization include the unit disc D with 
the Poincare metric, bounded symmetric domains, strictly pseudoconvex domains 
with metrics having reasonable boundary behaviour, or, provided one allows not 
only holomorphic functions but also sections of line bundles as elements of Hh, all 
compact Kahler manifolds whose Kahler form is integral. In all these cases, the 
choice of the spaces (ll.7[) which works are the weighted Bergman spaces Hh = 
L^ ol (fl, e - *' h u n ) (the subspaces of all holomorphic functions in L 2 (fl, e~®/ h uj n )), 
where n is the complex dimension of f2 and $ is a Kahler potential for ui (so, for 
ins tance for the unit disc H h = L 2 ol (D, ^(1 - |z| 2 )< 1 /^)-2 dz )y S ee [KS], [BMS] 
or |AE1] for the details and further discussion. 

Though this seems not to have been recorded explicitly in the literature, 
the whole formalism also extends seamlessly to spaces of vector-valued functions. 
In physical terms, this can be interpreted as accommodating the internal degrees of 
freedom of the quantized system. Namely, replacing the spaces Tin and L 2 (Sl, dfih) 
by the tensor products Hh <8 and L 2 (fl, dpkh) <8> C N (which can be viewed as 
spaces of C^-valued functions on ft), one can define in the same way the Toeplitz 
operators , where now the symbol <p can even be allowed to be a (N x iV)- 
mafra-valued function on fi. It is a simple matter to check, however, that this 
Toeplitz operator is just the N x N matrix of Toeplitz operators on Hh, 

and it immediately follows that (|1.5p remains in force in this vector- valued situation 
whenever it holds for the scalar-valued one. In particular, for scalar-valued func- 
tions <j> (i.e. <fi(z)jk = Sjk4>(z) for some (f> : Q — ► C), one recovers (| 1 . 5|) completely, 
with the same cochains Cj. 

In this paper, we work out a formalism, based upon certain spaces of matrix- 
valued functions, which could be looked upon, in appropriate cases, as a possible 
different approach to the quantization of the internal degrees of freedom of systems 
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whose kinematics is denned on complex phase spaces fi. Moreover, the more general 
setting adopted here, in that points in phase space are replaced by N x N matrices, 
could potentially be used to describe systems defined over non-commutative spaces. 

In more concrete terms, our spaces Sjh will be suitable subspaces (actually, rather 
small ones, in terms of codimension) of the spaces L 2 (fl,dfi h ) (£> C N of C^-valued 
functions on certain domains ft in C nxNxN associated to ft in a natural way. 
(Here, as before, n is the complex dimension of f2 and N is related to the number 
of internal degrees of freedom.) The Toeplitz operators are again defined by the 
formula (jl.ip . only with P replaced by the orthogonal projection onto Sjh, and 
the symbol <p can now be allowed to be a C ArxA '-valued function on i~2. Finally, 
there exists a canonical unitary isomorphism i : Sjh — * H-h <8> C N (with Tih as in 
(|1.7j) ). This means that the quantum system defined on $3h can be thought of as 
one possessing N internal degrees of freedom and moving on the phase space f2. 

The following facts then emerge from our analysis. 

(a) To any function <j) on fl one can associate, in a canonical way, a function <p 
on f2. (For reasons which will become apparent later, functions <p that arise 
in this way will be called spectral functions.) For any two functions (/>, ip 
of this form, the corresponding Toeplitz operators turn out to be unitarily 
equivalent via i to eg) / and Tffi ® /, respectively, acting on Hh <8> C N . 
Consequently, (jl.5| must hold (with the same cochains Cj), and, in this 
sense, our quantization contains the original scalar-valued Berezin- Toeplitz 
quantization, as well as its vector-valued analogue obtained by tensoring 
with (and using only scalar- valued symbols <p), mentioned above. 

(b) Let the unitary group U(N) of order TV act on C nxNxN by 

Z u = (U*Z 1 U,U*Z 2 U,...,U*Z n U) with Z= {Z u ...,Z n ) eC nxNxN . 

The domain $7 is invariant under this action, and functions <p satisfying 
4>(Z U ) = U*<p(Z)U MU S U{N) will be called V -invariant. All spectral 
functions are [/-invariant, but not vice versa. It is then the case that for any 
[/-invariant function 0, the Toeplitz operator is unitarily equivalent via 

l to the operator ® I on TLh <8> C^, where 7r^0 € C°°(f2) is a certain 
"average" of (p over the internal variables (reminiscent of "spin averaging" 
in quantum mechanical scattering theory). Not surprisingly, the operator Tth 
behaves nicely as h — > 0; owing to this, for any two U -invariant functions <fi, i/> 
one obtains a semiclassical expansion of the product T^T^ of the form 

oo 
r=0 

for some uniquely determined spectral functions v r . Thus, in this sense, 
the internal degrees of freedom disappear in the semiclassical limit, as they 
should. Using the isomorphism t and the facts mentioned in (a), this can 
also be recast into the language of the traditional vector- valued quantization 
discussed before; note, however, that now we are able to quantize not only 
the scalar- valued functions (which we have seen in (a) to correspond to the 
spectral functions on $7) , but a much wider class of observables correspond- 
ing to [/-invariant functions. 
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(c) Finally, for completely general functions 0, ij) G C°° (17), the semiclassical ex- 
pansion of the product TT'Tv' in the usual sense (i.e., in the sense of (II. 5|) ) 
does not exist. (There may be one, but the cochains Cj are then no longer 
uniquely determined unless one requires that their values always be spec- 
tral functions, and then they are no longer local (i.e. differential) operators, 
but rather involve some kind of averaging over a sort of t/(iV)-orbit of Z.) 
Consequently, such functions lead to quantum observables that have no clas- 
sical counterparts. The following situation is thus seen to emerge: while the 
Toeplitz operator corresponding a general function (p could be a legitimate 
quantum observable, only those functions which are [/-invariant, and conse- 
quently involve the internal degrees of freedom only in a "controlled" way, 
admit a semi-classical limit. In other words, only observables kinematically 
related to the phase space have semi-classical limits. (Note that even in the 
case where internal degrees of freedom are absent, i.e., N = 1, the model 
of quantum mechanics being used here is one where the wave functions are 
defined on phase space and not on configuration space.) 

The whole approach is applicable to any phase space J7 C C n admitting the 
ordinary (i.e. scalar-valued) Berezin- Toeplitz quantization. At the moment, we do 
not know how to extend it from domains in C™ to manifolds. 

For the simplest case of f2 = C, corresponding to a free particle on the real line, 
the results above have been obtained in AE2\ . For the reader's convenience, we re- 
view, in Section [5] below, the necessary material from that paper (without proofs), 
as well as from its precursor jAEGj , where spaces of matrix- valued functions of this 
type were first introduced. The quantization procedure is spelled out in Section [3] 
Hidden under surface in all these developments are also certain vector- and matrix- 
valued analogues of some reproducing kernels and coherent states; these in fact 
make sense in several more general situations as well (even though the quantization 
procedure may not lead to physically meaningful theories). We describe these in 
the last Section |4j 

A word of clarification is, perhaps, in order at this juncture. We are not sug- 
gesting here that the current formalism be used to replace the traditional quantum 
mechanical setup for describing systems with internal degrees of freedom. As far 
as traditional quantum mechanics is concerned, the present formalism, with wave 
functions described over matrix domains is an interesting alternative to it. Be- 
sides being well- adapted to studying the semi-classical limit, the present formalism 
can also be easily employed to build "quantum systems" which show no limiting 
semi-classical behaviour at all! One might venture a guess that such quantum sys- 
tems (which even admit a proper probability interpretation on "phase space" ) could 
point to some underlying non-commutative geometry. 



2. The case of the complex plane 

For the reader's convenience, wc briefly review here the salient facts from AE2 
and [AEGJ, which correspond to the simplest case of the quantization on fl = C. 

Consider the domain ft = {Z e C NxN : Z*Z = ZZ*} of all normal matrices 
in C JVxAr . By the spectral theorem, any Z £ Q can be written in the form 



(2.1) 



Z = U*DU, 
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with U G U(N) unitary and D diagonal; D is determined by Z uniquely up to 
permutation of the diagonal elements, and if the latter are all distinct and their 
order has been fixed in some way, then U is unique up to left multiplication by 
a diagonal matrix with unimodular elements. Consequently, there exists a unique 
measure dfi h (Z) on ft such that 

(2.2) [ f(Z)d» h (Z) = (Trh)- N [ f f(U*DU)e~^^ h dUdD V/, 
Jn Ju(N) Jc N 

where dU is the normalized Haar measure on U(N), dD is the Lebesgue measure 
on C , where we are identifying the diagonal matrix D — diag(di, . . . , d^) with 

the vector d={ d 1 ,.. .,d N )e C N , and ||D|| 2 = ||d|| 2 := |di| 2 + h \d N \ 2 . It can 

be shown [AEG] that 

(2.3) f Z* 3 Z k d f i h (Z) = S :jk k\h k I, 

Jn 

so that the elements 

(2.4) j = l,. ..,N, ^0,1,2,..., 

where xi, ■ ■ ■ , Xn is the standard basis of C N , are orthonormal in L 2 (tt, dfi h )®C N . 
Let S^h be the subspace spanned by these functions. 

In analogy with the scalar- valued situation, we next define for any <fi £ C°°(f2)<g) 
qNxn .j-kg Xoeplitz operator on $3h by the recipe 

T% ) f = P h {4>f), 

where Ph '■ L 2 (Q, dfi h ) — > Sjh is the orthogonal projection. Note that the last 
formula implies that \\T^ h) \\sj h ^sj h < \\<f>\\oo ■= sup Xen \\cj)(X)\\ C N^ C N . 

A C A ' xAr -valued function 4>(Z) of Z £ Q will be called spectral if it is a function 
of Z in the sense of the Spectral Theorem for matrices: that is, if there exists a 
function <j> : C — > C such that <p — 4^ where 

(2.5) <p*{Z) :=U- di&gj (<t>(dj ))•[/* if Z = U ■ diag^d,-) • U*. 

Further, as was already mentioned in the Introduction, the function will be called 
U -invariant if 

(2.6) <j>{u*zu) = u* 4>{z)u vu eu(N)yz en. 

Clearly, a spectral function is [/-invariant, but not vice versa: an example is the 
function tj>(Z) = |detZ| 2 /. 

The following results have been established in [AE2] . 

Proposition. ( |AE2j . Proposition 12) A function <j> is U -invariant if and only if 
there exists a function 4>{d\ \ cfej • ■ ■ , djv) from C x C w_1 into C, symmetric in the 
last N — 1 variables c^, • • • , djy, such that <f> = , where 

(2.7) 4>*(U*DU) := U* • dia gj (<t>(dj; d x , ....</, d N )) ■ U. 

The function <ft is uniquely determined by <fi. 

Further, (p is spectral if and only if <f> depends only on the first variable, i.e. if 
and only if </>(di; c?2, ■ ■ ■ , div) = 4>{di\ 0, . . . , 0). 

(Here we are using the notation <j># both in the sense of (|2.7[) and (|2.5[) . but 
there is no danger of confusion.) 
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Theorem. ( AE2], Theorem 10) If (j> = <j>* and tp = tp* are two smooth spec- 
tral functions on 17, then there exist unique spectral functions v r , r — 0, 1, 2, ... , 
such that 



in the sense of operator norms (i.e. as in (1 1.6} ). In fact, 

v r = C r (f,g)*, 

where 

(2.8) Cr (f,g) = -d r f -Wg 

r\ 



are the operators (| _?.#[) for n = 1 . 

Theorem. ( AE2J, Theorem 16) For a function f on C N and h > 0, let Khf be 
the function on C defined by 



TThf(zi) := / /(*!, z 2 , ...,z N ) e -(l«.l 3 +-+l«^l )A 



dz2 ■ ■ ■ dzjy 



Let <p = <P^ , "0 = ip be smooth U -invariant functions on f2 such that the partial 
derivatives of (f> and ip of all orders are bounded, and let C r be the bidifferential 
operators (| Then 



T (h)rp(h) ^ ST^ U rp(h) 



(2.9) 



r=0 

r U,) 

r=0 

in the sense of operator norms, where 



j',W>0, J ' ' ' 

where A' denotes the Laplacian with respect to the last N — 1 variables z%, . . . , zn, 
andf\z) :=/(z;0,...,0). 

Finally, for functions which are not [/-invariant, things seem to go wrong regard- 
ing quantization: namely, there is evidence that in general the semiclassical expan- 
sion of the form (|1.5p either does not exist, or if it exists then the cochains Cj have 
rather pathological properties (for instance, are not local operators — the value 
of Cj(4>,tp) at a point Z need not depend only on the jets of <p and tp at Z). 
In more detail: first of all, there exist functions <f> (even very nice and [/-invariant 
ones — for instance, <j>(Z) = \det{Z)\ 2 e~ Tr ^ z ' z h) for which \\T^ h) \\ -> as h \ 0; 
as a result, the cochains Cj in p.5[) have no chance of being uniquely determined, 
unless they are subjected to some additional condition. The only such condition 
which gives the right answer for spectral functions seems to be that Cj take values 
in spectral functions; let us therefore assume that this is the case. Second, there 
exist families of elements ' € $3h, labelled by Z G tt and x € (interpretable 
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as normalized reproducing kernels, or vector coherent states — see Section |4] below 
for more information), such that as h SO, there are asymptotic expansions 

oo 

r=0 

oo 

r=0 

for some C Arx ^-valued functions l r [<f>] and m r [</>, i/>] on SI, r = 0, 1, 2, If (|1.5[) 

holds, then we must therefore have mo[4>, if)] = Iq[Cq{4>, x/p)]. For spectral functions, 
lo turns out to be just the identity operator; since we have agreed that Cj take values 
in spectral functions, it follows that Co(</>, tp) = mo[<p, if}]. Now computations show 
that mo[4>, ip]{Z) is given by a rather complicated expression involving integration 
over the whole orbit {U*(dj(-, Xj)Xj)U} of the spectral projections dj (•, Xj)Xj> 3 = 
1, . . . , N, of Z under the unitary group U (N). The reader is referred to Sections 4-6 
of [AE2] for the full story. 

The appearance of $ and ifi , and not 4> and "0; m (|2.9[) means that the C" -1 part 
of <f) disappears in the semiclassical limit h — » 0, and only the projection (/^ , which 
lives on C, survives; that is, only the "spectral component" of the corresponding 
[/-invariant function <fi on SI. As mentioned before, all this means that we are 
dealing here with a quantum system which has N internal degrees of freedom, and 
that the full set of quantum observables of this system includes those which do not 
have classical counterparts, while even for those having the classical counterparts, 
the internal degrees of freedom — being purely quantum in this case — do not 
survive in the semi-classical limit. 



3. General domains 

We proceed to describe how the spaces from the preceding section can be adapted 
from the complex plane to any phase-space 17 C C™ admitting the ordinary (scalar- 
valued) Berezin-Toeplitz quantization. 

The appropriate matrix domain is 

Sl={Z = (Z 1 ,Z 2 ,...,Z n )EC n * N * N : 

ZjZ). = Z^Zj Vj, k = 1, . . . , n, and cr(Z) C fl}, 

i.e. the set of all commuting n-tuples of normal NxN matrices whose joint spectrum 
cr(Z) is contained in il. In other words, this means that in the decomposition (|2.1|) 
for the entries Zj, the unitary parts will be the same for all j: 

Zefi Z=(U*D^U,...,U*D^U), 

with U e U(N) and D {1 \ . . . , D (n) diagonal, 

and, if we denote the diagonal entries of the matrices by dJ (j — 1, . . . , n, 
/•• I VK 

(4 1) ,...,4 B) ) en > vfc = i,...,iv. 

With this notation, we define the measure on SI by 

1 N 
d ^( z ) = — row-i dU II d ^( d fe)> 
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where d k := (ujL , ■ ■ ■ ,d k ). We will also sometimes use the shorthand 

diag(di, . . . ,d„) 

to denote the n-tuple of diagonal matrices (D^ 1 ', . . . , D^) G f2. 

It remains to define the spaces jjj/j. Observe that since the Zj commute and 
cr(Z) G SI, the spectral theorem implies that, for any function / : VL — ► C, we can 
form the matrix f(Z u . . . , Z n ) -. f*(Z) G C NxN : specifically, 

(3.1) /#(Z) = C/*diag fe (/(4 1) ,...,4 n) ))C7 =U*di agk (f(d k ))U. 

Functions on SI of this form will be called spectral functions. We now define spaces 
Sj h C L 2 (ft,dfx h )<g>C N as 

S) h = span{/#(Z) X : / G Ll ol (Q,d» h ), X € C N }. 

Finally, recall also from the Introduction that a function <fi : fl — > C NxN is 
called U -invariant if ^(Z^) = U*d>(Z)U for all Z G tl and [7 G U(N), where 
Z C/ _ (u*ZtU, U*Z 2 U, U*Z n U). (Clearly, this reduces to the definition from 
Section [2] if n = 1.) 

Our main result is the following. 

Theorem 1. (i) The mapping 

(3.2) i:f(z)® X ^f*(Z) X 

is a unitary isomorphism o/£^ ol (f2, fih) ® C w onto Sjh- 

(ii) Under this isomorphism, the Toeplitz operator for a spectral func- 
tion <p = 4^ , corresponds to the tensor product ® I of the scalar Toeplitz 
operator T^' on L^ ol (f2, dfih) with the identity operator on C N . (In other words 

- to the Toeplitz operator T^j on L^ ol (fl, d^h) ® C N with scalar matrix-valued 
symbol <pl discussed in the second paragraph after $1.8$ in the Introduction.) 

(iii) Consequently, if <p = <fi# and if) = are two smooth spectral functions 
on O, then there exist unique spectral functions v r> r = 0, 1, 2, . . . , such that 

00 

r=0 

in the sense of operator norms (i.e. as in (j 1 . ) . In fact, 

v r = C r (f,g)*, 

where C r are the cochains \1.5\ from the ordinary (i.e. scalar-valued) Berezin- 
Toeplitz quantization on SI. 

(iv) A function <p : £1 — > Q NxN is U -invariant if and only if there exists a 
function 0(di;d2, . . . , djy) from SI x Sl^ 1 into C, symmetric in the last N — 1 
variables d 2 , . . . , d^r, such that <p = 4>^ > where 

(3.3) 0#(Z) := U*(dmg k (cp(d k - d 1; . . . ,d fe , . . .,d N )))U. 

The function 4> is uniquely determined by d>, and <p * s spectral if and only if <f> 
depends only on the first variable, i.e. if and only if (f>(di; d 2 , . . . , djy) = (f>(di). 
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(v) For a U -invariant function <p — 4># , the Toeplitz operator corresponds, 
under the isomorphism (jff.gj) . to the tensor product ® I, where 

1 f N 
-K h f{zi):= — , mjv _ 1 / f(z 1 ;z 2 ,...,z N )T[dn h (z 3 ). 

(vi) Consequently, for any two smooth U -invariant functions <p — 4>^ > V' = V 1 ^ 
on l~i smc/i i/iai the ordinary (scalar-valued) Berezin-Toeplitz quantization on fl is 
applicable to (f> and ip, 

oo 

(Q A\ rpWrpih) ^ , r rp(h) 

in the sense of operator norms, where C r are the cochains (| 1.5\ from the ordinary 
Berezin-Toeplitz quantization on Q. 

(vii) Finally, if, in addition, Q is one of the domains mentioned in the paraqraph 
after (| 1.8$ in the Introduction (examples of domains on which the scalar-valued 
Berezin-Toeplitz quantization is currently known to work), and 4> and ip have com- 
pact support, then \3.4-\ can be converted into an asymptotic expansion in pow- 
ers of h, i.e. there exist uniquely determined spectral functions v r , r = 0, 1, ... , 
such that 



(3.5) T^Tf^Y. hTT i* 



r=0 

in the sense of operator norms. 



The hypotheses in the part (vii) are made only for technical reasons, and can 
probably be weakened or dropped altogether. 

Proof, (i) For any /, g £ Ll ol (tt, d^ h ) and x, V € C N , we have 
(f*(Z) X ,g#(Z) V )= f V *g#(Zyf#(Z) X d» h (Z) 

= f i h (n) 1 ~ N f [ V *U*dmg k (gJ^)dmg k (f(d k ))UxdU]JdMd J ). 

JU{N) J Q N j 

Since, for any matrix X, 



(3.6) 

IU(N) 

we can continue the computation by 



[ U*XUdU=^I, 
Jmn) N 



N J ^ N k j 

1 ^(Q) 1 ^ rf X J2 ( / d ^ h ) N ^ f 9 Jd^f(d k )d fMl (d k ) 



^rf X N(f,g) 



= (x,v)(f,g), 

and the claim follows. 
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(ii) For f,g £ L^ ol (Q, dfih), XtV G and any function <j> on f2, we have by a 
similar computation as in (i), 

= ^ h (n) 1 - N f f 7 ? *C/*diag fc (^))diag0(d fc ) 
Ju(n) Jn N 

di^ifid^UxduHd^id,) 

j 

= ^ V*xN(<l>f,g) 
= {x,v){T^ h) f,9). 

(iii) follows immediately from (ii) and the ordinary Berezin-Toeplitz quantization 
on fi, upon tensoring with C N . 

(iv) Let : Q -> C A ' XJV be a [/-invariant function, and let D = (£>W, . . . , flW) 
be an element of f2 whose entries are diagonal matrices. For any complex numbers 
ei, . . . , ejy of modulus one, consider the matrix e = diag(e 1; . . . , ejy). Then e G 
U(N) and eZ?e* = Z? for any diagonal matrix D, whence D e = D; thus by the 
{/-invariance condition, 

0(D) =e* 0(D) e Ve Xj . . . , e N G T. 

Consequently, 0(D) is also a diagonal matrix. Define the functions /i,...,/jv : 
fi w -> C by 

/j(di;d 2 , . . -,d N ) := (0(D))^ where D = diag(di, . . . ,d N ). 

For any permutation a of the set {1, 2, . . . , TV}, let F a denote the permutation 
matrix [F a \jk — 5 a (j),k- Then F a G U(N) and 

F a T>F* = diag(d CT(1) , . . . ,d a{N) ) if D = diag(d x , . . .,djv). 

Thus by the [/-invariance condition again 

/<r(j)(di; d 2 , . . . ,djv) = /j (do-(i) ; do.( 2 ) > ■ ■ ■ ,d-<r(N))- 

It follows that fj is symmetric with respect to the last N — 1 variables di, . . . ,dj, 
...,d N and = /# for / = /i. 

Conversely, it is easily seen that any function of the form (|3.3|) is [/-invariant, 
and / # = g* f = g. 

Finally, the assertion concerning spectral functions is immediate upon comparing 
and (l3"TD . 

(v) Using (|3.3[) . the assertion (v) now follows by a similar computation as in the 
proofs of (i) and (ii): for f,g G L\ ol {Q,,d^hh) 1 x^V € and any : — * C 
as in (iv), we have 

(T$f#X,9*v) = J V*9 # (Zrd>*(Z)f#(Z) X dn h (Z) 



BEREZIN-TOEPLITZ QUANTIZATION 11 



= f x h (Cl) 1 - N / / r,*U*di agk (g(d k )) 
Ju(N) Jn N 



(diag fc (0(d fc ;d 1 , . . . ,d fc , . . . , d N ))) di&g k (f(d k ))UxdU JJtfyt/^d,,) 

3 

= 1 iihiSl) 1 -" V *X I Vip^^(d fc ;di > ... > d fcl ...,d JV )/(d fc )TTdMh(dj) 
= ti h (Q) 1 - N rf X I SpiMdi; d 2, • • • , djv)/(di) TT d^ h (dj) 

: / ff (di)( / #d i; d 2 ,...,d w ) TT^(d i ))/(d 1 )rf^(d 1 ) 



J=2 



= V * X f 5(dT)7r h ^(di)/(di)dA*A(di) 
Jo 

= V*x{(nh<f>)f,9) 

(vi) With (v) in hands, we obtain from the ordinary Bcrczin-Tocplitz quantiza- 
tion on f2, for any 0, t/j : — > C as in (iv), 

(h)rp(h) ^ rp(h) rp(h) j 

oo 

C r (Tr h rp,iT h if>) " 

r=0 

oo 

in the sense of operator norms (the last isomorphism being the one for spectral 
functions from part (i)), which proves (vi). 

(vii) Finally, to convert the last expansion into one of the form (|1.5p (i.e. in 
powers of h), we only need to exhibit a uniform asymptotic expansion for 7r^0, 
i.e. show that 

oo 

(3.7) i.^^^U 

r=0 

in the sense of norms in L°° (f2) , for some linear operators L r acting from functions 
on Q. N into functions on il. Indeed, since C T are bidifferential operators with smooth 
coefficients and <p,ip are assumed to have compact support, it will then follow that 

j,k>0 

in the sense of L°° norms, and in view of the inequality ||T« || < |M|oo, we can 
"apply T^ h > to both sides. 

In order to prove (|3.7p . it suffices in turn to show that there is an expansion of 
that form for 

(3.8) h k ^ h {yi) N ~ 1 u h (j){z) = h k I 4>(z- 1 z 2 ,...,z N )dfi h (z 2 ) ... dfi h (z n ), 
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for some k > 0, with leading coefficient which is positive on when <j) is identically f . 
Indeed, specializing this to <f> the constant one and dividing the two expansions 
gives (|3.7p . (The leading coefficient is needed to make sure that we are not dividing 
by zero.) 

Finally, for the situations where the ordinary Berezin-Toeplitz quantization is 
nowadays known to work (as summarized in the paragraph following {L8| in the 
Introduction), the measures [ih are taken to be e _ */' l d/i, where $ is a real-valued 
potential for the Kahler form u> and dfJ,(z) — uj(z) n — det[dd$(z)] dz —: g(z) dz is 
the Liouville measure. However, in that case the right-hand side of (13.81) reduces to 

/ fa; Z2, ... , Z N ) e -(*^) + -+*(^))/" g( Z2 ) . . . g( ZN ) d Z 2 . . . dz N , 

which has an asymptotic expansion of the desired form by the usual stationary 
phase method (or, rather, Laplace's method), with k — (N — l)(n— (the dimension 
of the variety on which — $ attains its global minimum)), and the leading coefficient 
being essentially the integral of <j>(z; Z2,---, ZN)g(z%) . . . 9{zn) over that variety with 
respect to the corresponding Hausdorff measure; see e.g. Fed , §4 of Chapter II, 
Hrm , Section 7.7, or MeJ, Chapter 7. This completes the proof. □ 

From a physical point of view, the sort of domains and Hilbcrt spaces envisaged 
in the above theorem could be used to describe system having n kinematic and N 
internal degrees of freedom. 



4. Some related reproducing kernels 

The original motivation that led the authors to the spaces like Sjh above did 
not actually come from quantization, but rather from an attempt to generalize 
to various vector- and matrix-valued setups the multifarious existing notions of 
coherent states from quantum optics (see e.g. |AAG] ). In our case here, these are 
given essentially by the "normalized" reproducing kernels of the respective spaces, 
see [AEG] ; for instance, as shown [AE2] . for the spaces S^h of Section [2] they are 
just the family k% ' , of elements of S)h, indexed by Z € ft (= the set of all normal 
N x N matrices) and vectors x G C , given by 

fcW (X) = K^{X,Z)K^{Z,Z)- l '\, 

where 

fc=0 

is the reproducing kernel of Sjh- 

Note that, despite the isomorphism from part (i) of Theorem [1] and the fact 
that the reproducing kernel for the corresponding space Hh is well known to be 
simply e( x ' v " h , the reproducing kernel (|4.I[) cannot be evaluated in a closed form, 
since the matrices X and Y* do not commute. For the same reason, it is impossible 
to evaluate in closed form any of the kernels from Theorem[T]even if the correspond- 
ing kernels for Tih are known. All one can do is to write them again in the form 
(|4.I|) , only the monomials need to be replaced by some general orthonormal basis of 
the space: namely, if {V'fc} is an arbitrary orthonormal basis of Tih = Z^ ol (fi, d/ih), 
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then by the well-known formula of Bergman Be the reproducing kernel of TLh is 
given by 

oo 
fe=0 

(the sum of the series does not depend on the choice of the orthonormal basis). 
In view of the isomorphism (|3.2p , it therefore transpires that 

oo 

(4.2) K™(X,Y) := £>#(X)^(Y)* 

k=0 

is the reproducing kernel of the space Sjhi in the sense that it has the reproducing 
property 

(4.3) /(Z) = f K (h \Z,X)f(X)d» h (X) VZ6fi,/6% 

JO 

In this short section we want to call attention to some situations when the quan- 
tization procedure from Theorem [1] does not apply, but there still exists a formula 
for the reproducing kernels like (|4.2|) . They all arise as Cartesian products of the 
spaces from Theorem Q] for the complex plane C, the unit disc D, and, more gener- 
ally, any one-dimensional domains £1 C C for which the ordinary Berezin-Toeplitz 
quantization works; for simplicity of ideas, we describe the space corresponding 
to C™ (the construction for the Cartesian product of any other n spaces of the 
above-mentioned type contains no additional new ideas). 

The domain in this case will consist of all (not just commuting) n-tuples of 
normal matrices: 

n = {Z = (Z x , ...,Z n )e C nxNxN : Z*Z 3 = ZjZ* Vj}. 

For the measure we take 

df*h( z ) = dti h {Zi) . . . dfi h (Z n ), 
where, abusing the notation a little, the d\i h on the right-hand side stand for the 
measure (|2.2j) on the set of all N x N normal matrices from Section [2] 

Finally, we define the space $3h to be the span in L 2 (Cl, dfi h ) (g> C N of the 
functions 

Z^Z^ . . . Z^Xj, fcl, • • • , K > 0, j = 1, . . . , N. 
Theorem 2. The mapping 



is a unitary isomorphism of fj^ onto L^ ol (C™, e _ " z " / h (nh)~ n dz\ . . . dz n 
Consequently, the reproducing kernel of S}h is given by 



K^ h \X,Y) 



E A 1 ■ ■ • A n Y n ■ ■ ■ Y l 
k 1 \...k n \h k ^+-+ k " 

fci,...,k„=0 



Proof. Let Zj = UfDjUj be the spectral decomposition (|2.1[) of Zj. Then by the 
definition of dfi hl 

(Z* Z k / . . . Zfrx, Z'l 1 Z* . . . Ztn) = j rfZ%« . . . Z*^ Z^ Z k 2 > . . . Z k ^ X d» h (Z) 
= (nh)- nN [ ... f [ ... [ f ] *U*D* n ^U n ...U^D* 1 jl U 1 

JC N JC N JU(N) JU(N) 
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UfDfrUi . . . U*D^UnxdUx ...dU n e --aimr>i+-+D*D n )/h dDi dDn 



Applying successively (|3.6[) to U%, U 2 , ■ ■ ■ , U n , we obtain 
(irh)- nN 

IC N JC N JU{N) JU(N) 



f ... / f ... f v*KK jn u n .-- u;ur u 2 Tt(d*^ d£ ) 

JC N JC N JU(N) JU(N) 



N 

UZD* 2 U 2 ... U*D k " U nX dU 2 ... dU„ e- ^(DtDi+-+D* n D n )/h dDl dDn 
{irh)- nN 



N 2 



{ ...j f ...f r,*U*D^U n ...UlD*^U 3 Tr(D* 1 j iD k 1 l ) 

JC N JC N JU(N) JU{N) 

Ty(D^D^) U£D k 3 s U 3 . . . U*D k *U nX dU 3 ...dU n 

e -W 1 D 1 + ...+D-D n )/h dDl ...dD n 



(iTh) 



N n 



f ... f r 1 *Tr(D^D k ^)...Tr(D^D k n ")x 
Jc N Jc N 

e -W 1 D 1 +...+D' n D n )/h dDl ...dD n 

v * x ff Mf^ f Ti{D*^D k -)e-^ D ' D ^ h dD 
"1 N JcN 
n 1 ,\-n N r 



m—l 

-N 

n x 

m=l ~ ' / j 
n -, AT 



1 



m=l 1=1 



m—l 

= ( x ,v) (z k >zt...z k %ztz?...zt). 
This settles the first claim. Besides, it shows that the function (X, Y) satisfies 

/ i^>(X, Z)Z*Zp . . . Z k " X d» h (Z) 

E_ klL X^:-^ J n z^...zrz k ^...z k ^ h{ z) 



jl,...,jn=0 



~ L> k\...k, ■ ' 1 •••*./•*•*•:••••'••<••" 

h,...,j n =0 
— Aj . . . A„ , 

i.e. has the reproducing property (|4.3p for functions of the form Z kl . . . Z kn x\ since 
the latter span all of Sjh by definition, the second part of the theorem also follows. 

□ 
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Although the reproducing kernels and the isomorphism (13. 2j) work out fine, what 
breaks down is that part (ii) of Theorem [TJ the Toeplitz operators on Sjh do not 
correspond, under the isomorphism above, to the Toeplitz operators on the Segal- 
Bargmann space L^ ol (C™, e~" z " / h (Trh)~ n dz\ . . . dz n ). The reason is the noncom- 
mutativity of Z\ , . . . , Z n — the reader can try to go through the beginning of the 
proof of part (ii) of Theorem [T] to see what is happening. For the very same rea- 
son, it also not possible to define spectral functions (and much less to describe the 
[/-invariant ones), and it is totally unclear at the moment how to achieve anything 
similar to the quantization from the previous section. 

We take this occasion to remark that one lands in even greater difficulties if 
one tries to deal with domains of arbitrary (rather than just normal) matrices. 
See Section 4 of ;AE2_ for details. 

On the other hand, it is clearly possible to use any other ordering of the entries 
of Z in Theorem [2] than Z\ . . . Z n (for instance, Z n . . . Z\)\ we omit the details. 
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